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ON  THE  STRUCTURE  OF  A LIAPUNOV  FUNCTIONAL 


FOR  A DIFFERENCE-DIFFERENTIAL  EQUATION 
WITH  ONE  DELAY 


Abstract:  A quadratic  positive  definite  Liapunov  functional 
that  yields  necessary  and  sufficient  conditions  for  the  asymp- 
totic stability  of  the  solutions  of  the  matrix  difference- 

differential  equation  Jx( tl/=  Ax(t)  + Bx(t-t)  is  constructed 

A*  l_ 

A , 

and  its  structure  is  analyzed.  The  method  of  construction  is 
a natural  generalization  of  the  same  problem  for  an  ordinary 
differential  equation,  and  this  relationship  is  emphasized. 


1.  Introduction . 

In  this  note  we  construct  a Liapunov  functional  for  the 
determination  of  the  asymptotic  behavior  of  the  solutions  of 
the  linear  autonomous  matrix  difference-differential  equation 
with  one  delay 


x(t)  = Ax ( t)  + Bx(t-x),  t > 0, 


where  x(t)  is  an  n-vector  function  of  time,  A and  B are  con 
stant  n x n matrices  and  T > 0. 

In  this  note,  we  generalize  results  presented  in  [7]  for 
a scalar  equation;  essential  to  this  generalization  was  the 

•i  StCtIM 

study,  presented  in  [2],  of  the  existence,  uniqueness  and  struc- 
ture of  the  solutions  of  a special  functional  differential  equa-  
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2. 


tion  which  is  intimately  connected  with  the  construction  of 
the  Liapunov  functional  presented  here. 

The  problem  of  constructing  Liapunov  functionals  for 
equations  of  this  type  has  been  previously  considered  by  Repin 
[8],  Datko  [3]  and  Hale  [5,6],  This  study  is  in  the  spirit  of 
the  previous  ones,  but  is  more  specific. 


2 • The  Dif ference-Dif  fercntial  Equation. 

Denote  by  ( [a,b]  ,_#n)  the  space  of  Lebesgue  square 

in tegrable ^functions  defined  on  [a,b]  with  values  in  i5Pn,  and  for 

a fixed  T > 0 consider  the  Hilbert  space  & = x L_  ( [-x , 0]  , 9Pn) 

T f°  T 2 

with  inner  product  cu^u^  = v^  + ^ (0 ) <P2  (0 ) dG  , where 

- T 

ui  = (v.,^)  £ i&,  and  the  naturally  induced  norm  is|  | ( v , <f> ) | ^ = 

T f ^ T n 

v v + <|>  (0)  (0)d0.  With  x:  [-T,00)  -*•  9t  , for  t > 0 we  denote 

— T 

by  xt  the  function  xfc:  [-i,0]  ->  9Pn , where  xfc(G)  = x(t+0). 
Consider  the  matrix  difference-differential  equation 


x(t)  = Ax (t)  + Bx(t-f),  t > 0, 


(2.1) 


where  A,B  are  n x n matrices,  x(t)  is  an  n-vector  and  t > 0, 
together  with  the  initial  conditions 


Xq  ( 0)  = T, , xQ  = <J>, 


(2.2) 


where  (£»<!>)  £ W. 


3. 


A solution  of  this  initial  value  problem  is  for  each  t > 0, 
a function  x G L2  ( [-T  , t]  , &n)  such  that  x is  absolutely 
continuous  for  t > 0,  satisfies  (2.1)  a.e.  on  [ 0 , t ] and 
x(0)  = £,  x(0)  = <j>(9)  a.e.  for  0 G[-t,0].  It  is  known  [1,6] 
that  (2.1)  — (2.2)  has  a unique  solution,  defined  on  [-t,»),  which 
depends  continuously  on  the  initial  data  in  the  norm  of  btf. 

The  initial  value  problem  (2.1)-(2.2)  can  be  rewritten  as 


d 

dt 


fxt(0) 


= srf 


xfc(0) 


(2.3) 


(xQ  (0)  ,xQ)  = (C,<f>)  € 


(2.4) 


where 


at 


xt(0) 

Axfc  ( 0 ) + Bx t ( ~ T ) 

ii 

3x  (0) 

Xt 

h — e ; °J 

(2.5) 


and  this  operator  has  a domain,  dense  in  b/f , defined  by 


SKjtf)  = {(€»*)  6 ^1  <J>  is  A.C.  in  [-t,0], 


4> ' 6 l2[-t,0]  , 4.(0)  = ?}• 
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The  operator  srf  is  the  generator  of  the  (^-semigroup  S/'it)  , where 
y{t)  : ^ •+  is  given  by  5^(t)  (£ , <|>)  = (x(t),xfc),  the 

solution  pair  of  (2.3),  (2.4). 

It  is  known  (1,6]  that  there  is  a constant  y such  that  the 
spectrum  of  srf  lies  in  the  left-half  plane  Re (A)  < y,  and  that 
for  every  E > 0 there  exists  a constant  K > 1 such  that 


I I I ( &,&■)  $ Ke 


(y+6)t 


(2.6) 


The  spectrum  of  srf  consists  of  those  complex  A which  satisfy  the 
characteristic  equation 


det[AI-A-Be"AT]  = 0. 


(2.7) 


Finally,  [5, 6],  a useful  representation  of  the  solutions  of  (2.1) 
is  given  for  every  t,  u > 0,  by 


[° 

xfc+u  (0)  = S(u)xfc(0)  + S (u-a-x)  Bxfc  (a)  da, 

J_T 


(2.8) 


where  the  matrix  S is  the  solution  of  the  matrix  initial  value 
problem 


ft  S(t) 


= S (t) A + S ( t-T ) B , 


= I,  S(t)  = 0 for  t < 0. 


(2.9) 


We  now  turn  to  the  construction  of  a Liapunov  functional. 
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3.  A Quadratic  Functional. 

Associated  with  the  functional  differential  equation 
(2.1)-(2.2),  or  (2.3)-(2.4),  and  motivated  by  the  results  in 
[7],  we  wish  to  consider  the  real  symmetric  quadratic  form  on 


V(£,4>)  = £Tm£  + e6x  f <J>T(0)Re26O<t>(O)de  + 

J -T 

+ CTQ(0)C  + 2CT  [ Q(a+T)e6(a+T)B<j)(a)da 


(3.1) 


+ 2 f ^)T(a)BTQ(3-a)e6  (a+6+2x)B^)(e)d6da  . 

~T  (X 


Here  6 is  a real  number,  M, R are  constant  n x n real  positive 

definite  matrices  and  Q(a)  is  a continuously  differentiable  matrix 
that  is  assumed  to  satisfy  the  initial  value  problem  for  the 
functional  differential  equation 

Q'(ct)  = ( AT+5 1 ) Q (a ) + e6xBTQT(x-a)  , 0 < a < t (3.2) 


Q (0)  = 0(0)-"  = Qq, 


(3.3) 


where  the  superscript  T denotes  transpose  and  is  a symmetric 

but  otherwise  arbitrary  matrix. 

Evaluation  of  the  Frechet  differentiable  Liapunov  functional 
(3.1)  along  the  solutions  of  (2. 3) -(2. 4)  with  initial  conditions 
in  ^(ja^)  yields  a function  of  time,  which  we  denote  by 
V(t)  = V (xfc (0) ,xfc) . This  function  of  time  is  differentiable  along 


% 


6. 

such  solutions,  and  a laborious  but  straightforward  computation, 
which  makes  use  of  (3.2)  and  (3.3),  shows  that  the  derivative  of 
this  function  along  these  solutions  is  given  by 

V(t)  = ^ V(xt(0),xt)  = -26V(xt(0) ,xt)  + (3.4) 

+ x£(0) [ (AT+6I)Q(0)  + Q (0) (A+6I)  + (AT+6I)M+  M(A+6I)  + 
+ e6TBTQT(T)  + e5lQ(x)B  + 2e6lR)xt(0) 

- e<5Th(xt  (0)  ,xt(x)  ) = U(xfc(0),xt), 

where 


h(xt(0) ,xt(-x))  = [x*(0) ,-e‘6Tx^(-T) ] 

A direct  application  of  Theorem  3.9  of  [9]  shows  that  for 
any  solution  with  initial  condition  on  {&,  V(t)  < U(xt(0),xfc). 

It  is  our  purpose  to  show  that,  through  the  functional 
(3.1)  and  its  derivative  (3.4)  it  is  possible  to  estimate  the 
rate  of  growth  or  decay  of  the  solutions  for  our  original  func- 
tional differential  equation  (2.3)  — (2.4) . For  this  purpose 
appropriate  choices  must  be  made  for  the  positive  definite  matrices 
M and  R,  of  the  constant  6 and  of  the  matrix  Q(a)  determined 
by  (3.2)  and  (3.3). 


R MB 


T 

B M R 


xt(0) 


-6t  , . 

-e  x (-x) 


.(3.5) 


L 


ftWiliiTiiiii  irivi1 1'  v 
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To  be  more  specific,  denoting  by  y = maxfRc  A|  det[Al-A- 
Be  ^T]  = 0),  we  wish  to  show  that  for  every  6>  0 and 
-6  = y + 2G  it  is  possible  to  choose  matrices  M,R  and  Q(a) 
satisfying  (3.2)  and  (3.3)  for  which  there  exist  positive  constants 
c^fC^  such  that 

clllU,<M||J,  < VU,$)  < c2|  | (£,*)  | d , (3.6) 


and 


V(Z,<P)  < -26V(5,4>). 


(3.7) 


These  last  relationships  imply  that  {V(5,<J>)  }l/'2  |'|  (C » ) fl  is  a 

' ' & 

norm  equivalent  to  the  original  norm  on  W and  that,  in  this  norm 


| (xfc(0)  ,xt)  | | < f|x0(0)  ,X(Jf 


& 


-6 1 
e 9 


(3.8) 


whereas  in  the  original  norm  on 


ll(xt<°).*t)lljr 


<V0,,V 


-6  t 


(3.9) 


These  estimates  are  precisely  those  stated  in  (2.6)  and  are  the 
best  possible.  It  should  be  noted  that  the  norm  induced  by  the 
square  root  of  the  Liapunov  functional  is  the  best  possible  one 
in  the  sense  that  it  yields  (2.6)  with  K = 1.  Moreover if 
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Y < 0,  then  the  Liapunov  functional  (3.1)  shows  that 
(2.1)- (2.2)  is  uniformly  exponentially  stable. 

First  of  all,  consider  an  appropriate  choice  for  the 

matrix  Q(a),  a solution  of  (3.2),  (3.3)  for  0 < a < x.  In  [2] 

it  was  shown  that  equation  (3.2)  with  initial  conditions  (3.3) 

has  a unique  solution;  moreover,  it  was  shown  that  the  linear 

2 

vector  space  of  all  solutions  of  (3.2)  has  dimension  n , and 

the  structure  of  these  solutions  as  well  as  simple  methods  of 
computation  which  take  advantage  of  the  particular  structure  of 
the  equation  were  presented. 

Associated  with  equation  (3.2) -(3. 3)  is  an  integral,  whose 
structure  is  motivated  by  a similar  integral  for  ordinary  differ- 
ential equation.  Indeed,  let  W be  a symmetric  matrix  and  let 
S (t)  be  the  solution  of  equation  (2.9,s.  Consider  the  expression 

#00 

Q (a)  = I ST(u)e6uWS(u-a)e6 (u_a,du.  (3.10) 

1 0 

Since,  for  every  6 > 0,  ||s(t)||  < Ke^Y+£^t  for  some  K > 1 and 
since  6 = -y  - 2G,  it  follows  that  this  integral  converges.  More- 
over, it  immediately  follows  from  the  definition  of  Q(a)  that 
Q(a)  = QT(-a),  Q(0)  = QT(0)  and,  since  S(t)  satisfies  (2.9), 
that  Q(a)  satisfies 





a 


-v,  ^ __  rp^.  -»  TiSTT  6 T 

Q • (0)  + q'T(0)  = (A  +61)  Q (0)  + Q ( 0 ) ( A+6 1 ) + B e Q (t)  + Q(x)e  B 

(3.12) 

= -ST(0)W  = -w. 


These  observations,  and  the  uniqueness  of  the  solutions  of 
(3.2)-(3.3),  show  that  Q(a),  defined  by  (3.10) , is  the  unique 
solution  of  (3.2)- (3. 3)  with  the  initial  conditions  prescribed  by 

the  second  equation  of  (3.11). 

It  is  easily  seen  that  the  map  W -*■  Q(0),  defined  by 


Q (0) 


S'1  (u)  e5  UWS  (u)  e^du, 

0 


as  a map  on  the  space  of  n x n symmetric  matrices  is  one-to-one, 
onto  and  it  maps  positive  definite  matrices  W into  positive 
definite  matrices  0(0) . 

With  this  particular  characterization  of  the  matrix  function 
Q(a)  it  is  possible  to  bring  into  evidence  the  particular  structure 
of  the  Liapunov  functional  (3. 1) .Indeed,  substitution  of  (3.10)  into 
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(3.1)  foe  Q(a)  yields,  after  some  rearrangements  and  interchange's 
of  integrals  , that  our  functional  is  of  the  form 


V(xt(0),xt)  = x£(0)llxt(0)  + e6 


T fe 


(0)Re2idx  (0)d 9 


+ j0Xt+u(0)Wo  xt+u(0)du’ 


(3.13.) 


Similarly,  using  this  notation,  equation  (3.4)  becomes 


U(xt(0),xfc)  = -26V(xfc (0) ,xfc)  + x£(0)[-W  + (AT+6I)M  + m(A+5I)  + 


ST  St  T -6t  T ^ ^ x^(0) 

+ 2e  R)x  (0)  -e  [xMO)  ,-e  TxM-t)] 

Z Z Z T -S  t 

B M R -e  xt(-T) 


. (3.14) 


Given  the  nonnegative  nature  of  the  last  term  in  (3.13)  it 

follows  that,  letting  A . (A)  denote  the  smallest  eigenvalue  of 

mm 

a symmetric  matrix  A, 

minU.  (M)  ,o~l6lTA  . (R) ) | | (€,♦)  I \h  < V(£,< p)  , 

nun  mm  ' ' - 


yielding  a value  of  c^,  for  equation  (3.6),  given  by 


c.  = mind  . (M),e  T)  in(Rl) 
1 mm  mm 


(3.15) 


11. 


Similarly,  since  from  (2.6), 


lXul0,ILn-  i Ke(Y+e)u||(x0(0),x0) 


ye 


& 


it  follows  that 


I r j X ( W ) p -y 

V(C,*)  < [max(  (M)  » e ° 1 T^max  tR)  1 + ^ 1 | | If. , »)  | ^ . 


This  yields  for  equation  (3.6)  a value  of  c ^ given  by 


I r [ X W n 

r,  6 t , , , max  „2 

c,  = max[X  (M) , e1  1 X (R) ) + ^ — K • 

2 max  max  26 


(3.1  0 


It  remains  to  be  shown  that  equation  (3.7)  holds.  For  this 

purpose  inspection  of  equation  (3.14)  indicates  that  it  is  necessary 
and  sufficient  to  show,  by  appropriate  choices  of  positive  definite 
matrices  W, R and  M,  to  have  the  matrix 


G = 


W - (A+6 1 ) Tm  - M(A+<5I)  - Re6T 


BTMe6T 


MBe 


6t 


Re 


6t 


(3.17) 


positive  semidef inite.  But  this  is  always  possible;  indeed,  a 
particularly  simple  choice  of  such  matrices  is  to  let  W,R  and 
M to  be  nonnegative  multiples  of  the  identity  matrix,  i.e. 

M = I,  R = kRI  and  W = kwI.  Letting  kR  = A (BBT)  and 


12. 


*w  = ■”^fO.»ma)((''+AT)  + 24  + 2o<ST/Amax(BUT)  ). 


then  (3.17)  is  positive  semidefinite  and  therefore,  for  these 
choices  of  M,R  and  W, 


V(xt(0),xfc)  < -26V(xfc(0) ,xfc)  (3.18) 

' 

for  every  (xt(0),xfc)  G 

It  is  to  be  noted  that,  with  these  choices,  the  Liapunov  func- 
tional is  particularly  simple  and  reduces  to 


260 

e dO  + 


V(xfc(0),xt)  = x£(0)xfc(0)  + e6T  /xmax(BBT)  1°  xj’(0)xt(0) 

— T 

+ max { 0 , (A+AT)  + 26  + 2e6Vx (BBT)  }.  (3. 19) 


1 ; 


max 


Xt+s(0)xt+s(0)e26Sds- 


max 


This  functional  is  a direct  generalization  of  that  used  in  (7),  and 
generalizes  those  used  by  Hale  [5,6]  and  Datko  [3].  We  recapitulate 
the  above  results  in  the  form  of  a 


. m*«*W~.*  ;^.r  , I 
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Theorem;  Consider  the  retarded  equation  x(t)  = Ax(t)  + Dx(t-i) 
and  the  Liapunov  functional  V given  by  equation  (3.1).  If 
Y — max { Re  A|  det[AI-A-Be  = 0}  and  G > 0,  then  there 

exist  constant  positive  definite  matrices  M and  R and 

T 

a differentiable  matrix  Q(a),  0 < a < t with  Q(0)  = Q(0) 
such  that  the  functional  V is  positive  definite,  bounded 
above,  and  V < 2(y+£)V. 


Of  course,  if  y < 0,  the  above  result  implies  exponential 
asymptotic  stability;  moreover,  the  rate  of  decay  is  precisely  the 
expected  one. 

4 . A Comparison  with  Ordinary  Differential  Equations. 

In  this  section,  it  is  our  object  to  point  out  the  intimate 
relationship  between  the  results  obtained  for  our  functional  equation 
and  the  classical  results  on  the  construction  of  Liapunov  functions 
for  ordinary  differential  equations. 

Recall,  [4],  that  given  the  system  of  ordinary  differential 
equations 


x(t)  = Cx (t) , 


(4.1) 


where  C is  an  n x n matrix,  a Liapunov  function  for  this 


14. 


system  can  be  always  taken  as  the  quadratic  form 

V (x)  = xPPx,  (4.2) 

where  F is  a positive  definite  matrix.  Moreover,  if 
y = max {Re  X|  det[Al-C]  = 0}  and  if,  for  G > 0,  -6  = y + 2G, 
then  given  an  arbitrary  positive  definite  matrix  W,  the  algebraic 
equation 

(C+6 1 ) TP  + P(C+SI)  = -W,  (4.3) 

has  a unique  solution  P which  is  positive  definite.  This  matrix  p, 
if  used  in  (4.2)  along  the  solutions  of  the  differential  equation 
(4.1)  yields,  upon  differentiation, 

<f(x(t))  = -2SV(x  (t ) ) - xTWx  < -26V  (x  (t)  ) . (4.4) 

Furthermore,  the  unique  positive  definite  solution  of  (4.3)  can 
be  obtained  as  the  integral 

P = f”BCTue6uWecueiudu.  (4.5) 

1 0 


Let  us  now  bring  into  evidence  the  relationship  between  the 
results  obtained  in  the  previous  section  and  these  results.  Th-_ 
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functional  differential  equation  under  consideration  is 


x(t)  = Ax  (t)  + Bx(t-x),  t > 0 . 


(4.1' ) 


Once  again,  we  assume  that  y = max{Re  \\  det[AI-A-Be  = 0} 

and,  for  G > 0,  let  -5  = y + 2G.  The  Liapunov  functional  is 
then  of  the  form 

V(xt(0),xfc)  = x£(0)Mxt(0)  + e*51  j x^  ( 0 ) Re2  (5°xt  ( 0 ) dO  (4.2') 

+ x£(0)Q(0)xt (0)  + 2x£(0)  Q(a+T)e6 (a+T)Bxt(T)dO  + 

-T 

+ 2 °x^(a)BTQ(3-a)e6(a+6+2T)Bxt(3)d3da. 

•'-T  a 


The  choice  of  the  positive  definite  symmetric  matrices  M and  R 
in  this  expression  is  rather  arbitrary.  Their  purpose  is  to  insure 
the  strict  positive  definiteness  of  the  functional  on  the  Hilbert 
space  i&;  it  should  be  noted  that  if  M = R = 0,  the  functional 
(4.2')  is  positive,  but  does  not  satisfy  a relationship  of  the  type 
V(xt(0),xt)  > c1|  | (xt(0) ,xfc)  | for  ^ > 0.  The  requirement  that 

the  matrix  G of  equation  (3.17)  be  positive  semidefinite  is  always 
satisfied  for  M = R = 0 and  W positive  semidefinite.  Given  an 
arbitrary  positive  definite  matrix  W,  it  is  always  possible  to 
select  positive  definite  matrices  M and  R so  that  G is  positive 
semidefinite. 


The  choice  of  the  continuously  differentiable  matrix  Q(a), 


0 < a < T is  critical.  Given  an  arbitrary  positive  definite  matrix 
W it  must  satisfy  the  functional  equation 


Q'(a)  = (AT+6l)Q(a)  + BTe6xQT (x-a ) , 0 < a < t, 

Q (0)  = QT(0), 


(4.3*  ) 


with  the  condition 


(AT+6 I ) Q ( 0)  + Q(0)(A+6I)  + BTe6TQT(x)  + Q(x)e6TB  = -W.  (4.3") 


In  the  previous  section  it  was  shown  that  such  a Q(a)  always 
exists  and  is  unique. 

With  such  a choice  of  Q(ot),  one  then  obtains  that,  along  all 
solutions  of  the  functional  differential  equation  (4.1') 


V (xfc (0) f xfc)  <-26V(xfc (0) ,xfc) . 


Moreover,  such  a matrix  Q(a)  exists,  is  unique 
of  it  is  given  by  the  integral 


Q (a)  = f ST(u)e6uWS (u-a)e6 (u_a)du, 
1 0 


where  S(t)  is  the  solution  of  equation  (2.9) 


(4.4') 

and  a representation 

(4.5') 
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The  strong  relationship  between  the  unprimed  and  primed 

equations  is  now  clear.  Indeed,  note  that  for  t = 0,  the  matrix 

C in  (4.1)  becomes  A + B and  the  matrix  P in  (4.3)  becomes 

Q(0)  + M and  all  the  primed  equations  become  the  unprimed  ones. 

Equations  (4.3')  and  (4.3")  are  of  a much  more  complex 

nature  than  the  familiar  algebraic  equation  (4.3).  However,  in 

spite  of  its  appearance,  the  linear  vector  space  of  the  solutions 

of  (4.3')  is  not  infinite  dimensional  but,  as  v/as  pointed  out  in 

2 

the  previous  section,  has  dimension  n . Hence,  although  the 
problem  of  construction  of  the  Liapunov  functional  for  the  func- 
tional differential  equation  does  not  reduce  to  the  solution  of 
an  algebraic  equation  such  as  (4.3),  it  reduces  to  the  solution  of 

a linear  differential  equation,  equation  (3.10),  which  can  be 
accomplished  without  difficulty. 
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